Magnetic Anisotropy in the Molecular Complex V15 



N.P. Konstantinidis 

Department of Physics and Department of Mathematics, Trinity College, Dublin 2, Ireland 

CN ■ D. CofJey* 

I Department of Physics, State University of New York at Buffalo, Amherst, NY 14-260 

^ ! (February 1, 2008) 

^ ■ We apply degenerate perturbation theory to investigate the effects of magnetic anisotropy in 

Q ' the magnetic molecule V15. Magnetic anisotropy is introduced via Dzyaloshinskii-Moriya (DM) 

, interaction in the full Hilbert space of the system. Our model provides an explanation for the 

' rounding of transitions in the magnetization as a function of applied field at low temperature, from 

\ which an estimate for the DM interaction is found. We find that the calculated energy differences 

of the lowest energy states are consistent with the available data. Our model also offers a novel 
explanation for the hysteretic nature of the time-dependent magnetization data. 
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PACS numbers: 75.10.Jm Quantized Spin Models, 75.50.Ee Antiferromagnetics, 75.50.Xx Molecular 
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I. INTRODUCTION 



' Magnetic molecules have attracted significant attention in recent years as systems where macroscopic quantum 
phenomena are displayed, such as relaxation by tunneling of magnetization through a potential well [1]- [4]. These 
molecules form crystals with very weak intermolecular interactions, thus experiments reflect directly the properties of 
'"^ individual molecules, which can be described by magnetic Hamiltonians. The most studied molecules are a cluster of 
■ manganese ions known as Mni2 [5], and clusters with iron ions, Fe„ [6-8]. 
O \ Mni2 and Fes possess strong easy axis anisotropy and their states can be characterized to a very good approximation 
I by < S*^ >, the expectation value of the total spin along the easy axis, z. Their ground state has a total spin S = 10, 

and to a good approximation they are described by a single-spin model, where only the S* = 10 multiplet states are 
. considered [9,10]. Restriction to the 5 = 10 multiplet reduces drastically the Hilbert space of the problem. Iron 
^ ' atoms also form ferric wheel clusters with predominantly antiferromagnetic interactions. In these the total spin of the 
\ ground state is 5* = [11]- [15]. Spin dynamics of Peg and Peg have recently been calculated by Honecker et al. [16] 
using exact diagonalization. 

Another molecule of the same class is a cluster with 15 Vanadium ions, known as V15 [17], each of which has 
5 = i. This is similar to some of the Pe„ systems in that the principal interactions are antiferromagnetic. This 
fS| system has also been investigated using effective spin models in which sets of spins are replaced by one effective spin. 
. However, effective three-spin models which contain only isotropic terms [18]- [20] can not explain the broadening of 
' the spin transitions as a function of magnetic field in equilibrium magnetization measurements, as was pointed out 
^ by Chiorescu et al. [18]. Miyashita and Nagaosa [21] considered a three-spin model for V15 which included a term, 
'^ij^i^j- This term is responsible for the quantum tunneling of magnetization states and for the smoothing out 
of transitions between magnetization states in applied magnetic fields. However, a more microscopic model is needed 
if a quantitative comparison with data is to be made. In particular, the tunnel splittings, which are important for 
Q magnetic relaxation, are very sensitive to Hilbert space truncation and the full Hilbert space of the molecule has to 
Q . be taken into account [22]. Here we solve for the lowest energy states using degenerate perturbation theory [23,24] in 
the full Hilbert space, which has 2^^ = 32, 768 states. 

Apart from the opportunity to study the magnetic properties of finite systems which these systems provide, Leuen- 
berger and Loss [25] have proposed the use of Mni2 as a qubit, a quantum bit, in a quantum computer. In this type of 
^ [ application it is important that prepared spin states do not decohere on the time scale of the operations on the system 
- - < necessary to read, write and carry out the manipulations necessary for a computation. Consequently it is important 
to understand the mechanisms for magnetic relaxation in these molecules. These mechanisms can be probed by the 
application of time dependent magnetic fields [26,27]. 

The magnetic hysteresis studies for V15 show so-called "butterfly" curves which are different from the ones of 
Mni2 and Peg [18,28]. This has been explained in terms of spin-phonon interaction using the Landau-Zcncr model for 
transitions between two states [29]. However, the low number of phonons at very low temperatures points to a different 
mechanism for magnetization relaxation, intrinsic to the molecule. Here we examine the role of magnetic anisotropy 
described by the Dzyaloshinskii-Moriya (DM) interactions as a relaxation mechanism. This interaction arises from 



1 



spin-orbit coupling and was introduced by Dzyaloshinskii phcnomcnologically by using symmetry considerations [30] . 
Moriya derived it microscopically by extending the Anderson theory of superexchange to include spin-orbit coupling 
[31]. The DM term arises in low symmetry crystals. The fact that finite systems have surfaces where symmetry is 
reduced ensures that the DM term is an important part of any microscopic description of magnetic molecules [32] . 
There are other anisotropy sources such as dipolar and hyperfine fields which have very small magnitudes, about 1 
mT and 40 mT respectively [18]. Magnetic anisotropy has also been modeled in Ising Hamiltonians by and S~ 
terms to reproduce the experimental data [9,10]. Single site anisotropy terms do not contribute in V15 since af = 1, 
i = x,y,z for S = ^. Therefore the DM interaction is the leading order anisotropy term for this system. 

The model Hamiltonian, which includes a Heisenberg exchange interaction and a DM anisotropy term, is solved here 
with degenerate perturbation theory [23,24]. It is found that inclusion of the DM terms quantitatively explains the 
widening of the spin transition in the equilibrium magnetization curve. The DM terms are determined quantitatively 
by comparing with experimental data of Chiorescu et al. [18], and the resulting splitting of the low energy groups 
of states is in agreement with experimental estimates, ^ 10^^ K. Furthermore, it is shown with a model calculation 
that the non-equilibrium magnetization curve is intrinsically hysteretic due to the presence of the DM interaction. 
This interaction leads to rapid oscillations wth frequencies given by its magnitude. The form of the hysteresis curve 
depends on when during these oscillations the magnetic field is reversed. 

The paper is organized as follows. In section II we discuss the magnetic description of the V15 molecule and in 
section III we introduce the model including the DM interaction determined by symmetry. In section IV we present 
the properties of the model and compare with experimental data. Finally in section V we give our conclusions. 

II. Vis 

The chemical structure of V15 is given by the formula Kq [V(5'As6042(H20)] • 8H2O [17,33,34]. There are fifteen 
V^^ ions with spin S = ^ placed in a quasi-spherical layered structure (figure 1). These spins sit on two hexagons with 
a triangle sandwiched between them and all interactions are antiferromagnetic in the simplest model [17]. The spins 
in each hexagon are paired via a strong coupling J = 800 K [34]. These pairs are connected with a nearest neighbor 
interaction J' = 150 K. There is also a diagonal interaction J" = 300 K. The triangle spins do not interact directly 
with each other. They interact with spins from the top and bottom hexagons, and the interactions are Ji = J' and 
J2 = J" . The molecule has trigonal symmetry and the space group is R-^^. The unit cell includes two molecules but 
intermolecular interactions are negligible. 

Since all interactions are antiferromagnetic the total spin of the ground state is S = \. This is in contrast with 
Mni2 and Fcg, where S = 10. In the absence of anisotropy the ground state is four-fold degenerate, and it is made 
up of two doublets. The lowest lying excited states form a quartet which is separated from the ground state doublets 
by « 3.7 K. The splitting between the lowest energy multiplets from the rest of the spectrum has been found from 
EPR spectra to be ~ 500 K [35]. Magnetic anisotropy leads to tunnel splittings between the two doublets, which are 
^ 10~^ K [36], in contrast with the very small values for the other two molecules 10~^ — 10""'^^ K). There is no 
macroscopic quantum tunneling for V15 according to the Kramers theorem [37] , since the system has an odd number 
of 5 = I spins and the total spin is half- integer [38], giving a doubly degenerate ground state. 

From equilibrium magnetization measurements, it is clear that the ground state has 5 = ^ up to a critical value 
of the field, where a transition to 5 = | is detected [18]. This is true for temperatures smaller than 0.9 K. The 
broadening of the spin transitions can not be attributed to temperature alone, and various factors like dipolar or 
nuclear hyperfine field distributions are too small, about 1 mT and 40 mT respectively [18]. This broadening can 
be explained by anisotropic DM interactions, which were also found to be important for the description of neutron 
scattering data and EPR measurements in Mni2 [39]. More recently, analysis of neutron scattering data on V15 by 
Chaboussant et al. [40] has shown the need to include DM interactions. Garanin and Chudnovsky [41] have suggested 
that the width of EPR lines could also be explained by dislocations in the crystal lattice. Here we assume that there 
are no defects in the crystals so that data reflects the properties of the V15 clusters. 

III. MODEL 

A microscopic treatment of the problem of the origin of magnetic anisotropy is considered using the full Hilbert 
space with 2^^ = 32, 768 states. This approach based on degenerate perturbation theory makes it possible to include 
the effect of spin-orbit coupling through the DM interaction. Initially perturbation theory is applied in the absence 
of magnetic anisotropy. The unperturbed Hamiltonian, Hq, takes into account the spin singlet-triplet correlations 
present due to the strength of the bond between alternating pairs of spins in the hexagons being more than twice as 
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large as any other coupling. The ground state of Hq is eight-fold degenerate due to uncoupled spins in the triangle 
(figure 1). Then the perturbation includes the remaining terms for bonds on which exchange constants are smaller. 
An effective Hamiltonian, H'^^^ , is generated for the degenerate subspace. Since singlet-triplet correlations of the 
prevalent J bond are built into Hq, the series expansions generated by perturbation theory are absolutely convergent. 
Diagonalizing H^f-^ gives the eight lowest energy states of the system which are separated from the next set of states 
roughly by the energy splitting between singlet and triplet states in Hq. The Hamiltonian is of the form 

H = Ho + XHi (1) 

where A is the perturbation parameter with 

Hq = J{Si ■ S2 + S3 ■ S4 + S5 ■ Sq + Sio ■ Sii + S12 ■ Si3 + Si4 ■ S15) (2) 

and 

Hi = J'{S2 ■ Ss + S4 ■ S5 + Si ■ Se + Sii ■ S12 + Sis ■ S14 + Sio ■ S15) + 

J" {Si ■ S3 -\- S'i • Sc, -\- Si ■ 5*5 + 5*10 • S12 + S12 ■ Sii + 'S'lo • 'S'14) + 

Jl{S2 • 5*7 + 5*7 • Sii + 5*4 • 5*8 + 5*8 • 5*15 + 5*6 • 5*9 + 5*9 • S13) -\- 

J2{Si ■ Si + Si ■ Sici + S3 ■ Ss -\- Ss ■ Si4 + 5*5 • 5*9 + 5*9 • S12) (3) 

The ground state of the unperturbed Hamiltonian is composed of a singlet at each of the six hexagon pairs. When 
A ^ the perturbation mixes in states where one or more of the pairs are excited and generate the elements of H"^-^-^ . 
Since J is a very strong bond, this lowest multiplet of 8 states will be well-separated in energy from the rest of the 
spectrum, and the low temperature behavior is well described by the eigenfunctions of H*^^^ . 

In equation (1) the exchange constants given in [34] are varied when A = 1 so that the gap from the doublets to 
the quartet agrees with the experimental value of « 3.7 K. The following parameters are kept: 

J=l, Ji = J' = ^J, J2 = J" = ^J (4) 
' 800 800 ^ ' 

The unit for J is 800 K and the energy gap to the excited state with these values is equal to Ai = 3.61 K. Thus, the 
low energy diagram of V15 is as shown in figure 2. Aq ~ 10~^ K is the tunnel splitting between the two doublets. This 
splitting is not described by the present Hamiltonian. A test for the success of perturbation theory is the calculation 
of the expectation values of the square of the total spin S"^ and its projection on the z axis, S^, for the calculated 
wavefunctions. These expectation values are half- integer for S^ and integer for 5^ at A = 1 with a precision equal 
to the one used to generate the perturbation expansions (double precision-15 digits). This confirms the success of 
perturbation theory. 



A. Magnetic Anisotropy 
The anisotropic interaction between two spins Si and Sj due to spin-orbit coupling is [30,31]: 

' Si X Sj + Si ■ F; j • Sj (5) 

The first term is the DM interaction and it is first order in spin-orbit coupling and antisymmetric. We refer to Dij 
as the DM vector. The second term is second order in spin-orbit coupling and symmetric, with Fij a second-rank 
tensor. Thus the DM vector provides the lowest order anisotropy term. 

As a minimal model, the DM interaction will be considered on the J bonds only, the strongest exchange constants 
in the system. The DM interaction is: 

Dij .SiXS,= Dl^iSfS] - S^S]) + Dl^iS^S] - S^S]) + DtJSfSy - SfS^) (6) 

Symmetry operations of the V15 molecule, rotations of ^ and ^ around the axis that passes through the centers of 

the hexagons {z), have to leave the form of the Hamiltonian invariant. This constraints the Dij to be related to one 
another. If 

Di,2 = X + y + z (7) 



3 



then these three parameters determine the DM interaction in the other two J bonds of the bottom hexagon as: 

^3,4 = \ i-DU + Dl,) X - i(V3 Dl, + Dl^) y + z 

D5,6 = -l {Dl^ + V3 Dl^) x+\{V^ Dl^ - Dl,) y + D^^ z (8) 

The sites 1, • • • , 6 arc those in the bottom hexagon. Due to absence of reflection symmetry with respect to the triangle 
plane the DM terms in the upper hexagon will differ from the ones in the lower one. However, to minimize the number 
of free parameters in the model, -Dio^i = ^i,2- Then from symmetry D\2,\z — ^5,6 and i5i4,i5 — -D3 4. 

In the same manner symmetry determines the elements of j . Considering the symmetric part of the spin-orbit 
coupling again only on the strong bonds (J), and if 

12 12 12 
pa:y Y-VV -pv'z 
J- 1,2 J- 1,2 i 1,2 
-pxz rizz 

J- 1,2 J- 1,2 ^ 1,2 

then symmetry determines Fij by: 

ri,j = u Fi,2 

where U is a. rotation by (f) around the z axis: 

(cost/) sin(j) \ 
—siruj) cos<p 
1) 

F3,4 and T^ q are determined in this manner and, as in the case of the DM terms, to minimize the number of free 
parameters in the model Fio,ii = Fi_2- Then from symmetry ri2,i3 = Fs^e and Fi4_i5 = F3^4. 
The DM term added to the Hi term in equation (1) is of the form: 

Hdm = ^1,2 -Six §2+ ^3,4 • S3 x ^4 + ^5,6 -55x56 + 

-Dl0,ll • •S'lO X ^ii + r>12,13 • 512 X 5i3 + r>14,15 ' 5i4 X 5i5 (12) 

This term is considered as part of the perturbation, Hi, and is scaled with the perturbation strength A. Prom now 

on = Df 2' Dy = D\ 2 and = Df 2- This minimum model is investigated for magnetic anisotropy and the 
parameters, D'' , and , are determined by comparing with experiment. 

We have found that the symmetric Fi j terms, which are second order in spin-orbit coupling, have little effect in 
the magnetization curves. Once we have determined the parameters of the minimal model discussed above, we will 
return to this point. In the case where the moments lie on an open linear array it is possible to introduce local axes 
so that the DM term, Dij ■ Si x Sj, is absorbed into an effective symmetric term [42]. However, the connectivity of 
the bonds between the spins on the V sites ensures that this term can not be transformed away in this manner. 

IV. RESULTS 

The effect of the DM terms in the energy spectra is investigated by varying the strength of D^, and D^. The 
average energies of the two doublets and the quartet are plotted for different values of the D vector in figure 3. The 
changes in the energy become more pronounced as |5| increases. However, one sees that even DM terms ~ 100 K lead 
to energy shifts ~ 20 K. The deviations from the average energies of the two doublets are plotted in figure 4, while 
for the quartet in figure 5. It can be seen that even when \D\ ~ 100 K the sizes of the energy splittings in the two 
doublets and the quartet are modest. This suggests that the energy level differences estimated from the experiment 
are comparatively insensitive to D so that it can not alone determine D. 

These spectra arise from the form of the effective Hamiltonian for the triange spins, H*^^^ , calculated with per- 
turbation theory. The effective Hamiltonian for the triangle spins is of the form H^ff = iJ*"^ + iJ-o//-*«9, where 
jjdzag jg constant diagonal matrix in the basis of total 5^ states which is three orders of magnitude larger than the 
terms in jjoff-d^ag ^^^^^ depends on the magnitude of the anisotropy term, as can be seen in figure 3. All the total 
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(9) 
(10) 

(11) 



S^- sectors of the triangle spins are coupled in H"^ f^''-'^"-^ even when the anisotropy term is ■ Si x Sj on the 
800 K hexagonal bonds, case (I) in figure 4 and figure 5. 

By contrast, when one considers a three-spin model with a Heisenberg coupling, Jq Si ■ Sj, and an anisotropy 
term, ■ '^■j Si x Sj, between the triangle spins, the quartet remains degenerate while the doublets are split into 
time-reversed pairs separated by \/3D^ . Unlike the result of this three-spin model, the energies of all quartet levels are 
changed by the anisotropy term, as can be seen from figures 4 and 5, and the splitting of the higher energy multiplet 
is larger than the low energy one. This points to the difficulty in understanding the effect of the anisotropy terms on 
H'^-l'^ in terms of a simple model for the triangle spins. More generally it supports the conclusions of De Raedt et al. 
[22] on the dangers of Hilbert space truncation. 

A. Magnetization 

The magnetization of V15 was calculated for various magnetic fields and temperatures as a function of the DM 
terms. In the absence of DM terms and at zero temperature, S^ is a good quantum number and the magnetization 
stays constant cixcept at critical fields where there are level crossings. In the first level crossing the ground state 
selects the state of the zero field ground state doublet with the higher expectation value of the spin, while in the 
second it switches from this doublet to one of the quartet states. When the temperature differs from zero, the jumps 
are smeared out and at high enough temperatures there are no sudden changes in the magnetization. 

When the DM terms are non-zero at zero temperature S"^ does not commute with the Hamiltonian, and the same 
is true for S^ if either or is non-zero. The magnetization along the field was calculated for = 50 K and 
— 100 K at a very low temperature, T = 0.001 K, in the case of a field along the x axis and a field along the z 
axis and is plotted in figure 6. The gyromagnetic ratio used is g = 2.0023 [43]. It observed that the magnetization 
depends on the degree of anisotropy, D^. In addition, it is seen that the jump is more smeared out as the magnitude 
of increases from to 100 K. Magnetic anisotropy is demonstrated by the dependence of the magnetization along 
the direction of the field on the direction of the field itself. 

The broadening of the spin transition can also be observed in figure 7. There the magnetization < S'^ > is plotted 
as a function of the magnetic field for four different temperatures and a DM term of magnitude \/3 100 K. The 
magnetization curve is almost identical for T = 0.001 K and T = 0.1 K. This shows that the broadening of the 
transition is not a temperature effect, but it is associated with the DM terms. For high enough temperatures the 
susceptibility is almost constant with the field. 

In figure 8 the calculated magnetization is plotted for three different values of with = = 0. It is seen 
that the width of the transition gets bigger with the increase of the DM terms. The origin of the width can be seen 
in the energy spectrum plotted in figure 9. The presence of matrix elements between the states in the vicinity of the 
transition field (~ 3 T) leads to level anticrossings. Pairs of degenerate states from the lower multiplet in zero field 
split up with one hybridizing with a member of the quartet multiple with i=a |, while the other has no matrix 
element with this state. 

The width of the magnetization transition is an indication of the effect of the DM terms, as can be seen from figures 
7 and 8. The transition is characterized by its midpoint BtransiUon and width- difference AB in magnetic field values 
where it has changed by 10 % from its low and high field average spin values on either side of the transition. The 
results for various DM terms are shown in table I, where it can be seen that BtransiUon depends sensitively on the 
form of the DM terms. This dependence on the DM terms is illustrated in figure 10 where we plot BtransiUon versus 
AB for two different forms of the DM vector, one in which ^ K, = = and the other in which 7^ K, 
]jx _ jjy _ Q rpj^g values of BtransiUon and the Corresponding values of AB shown are for different magnitudes of 
and at two different temperatures. Although the magnetization curves for these two temperatures appear to lie 
on top of one another in figure 7, there is a difference which should be taken into account in quantitative comparison 
with data. The effect of increasing temperature is seen to push BtransiUon to higher fields, while AB is seen to be 
non-monotonic. For small magnitudes of the DM term AB is larger at T = 0.1 K than at T = 10"'^ K but this is 
reversed for a sufficiently large magnitude depending on the two forms considered. 

These temperature dependences can be understood as follows. BtransiUon moves to higher applied fields with 
increasing temperature because the higher energy low magnetization state contributes significantly to the thermo- 
dynamic average for larger energy differences between the ground and excited state, that is for larger applied fields 
beyond the field at which the two states are degenerate. The temperature dependence of AB is more complicated. 
Both the temperature and the Dij vectors contribute to AB. Because the temperatures considered here are so low 
compared to the variation in the magnitudes of the Dij vectors, the difference in the smearing effects of the different 
temperatures, 0.1 K, can be neglected compared to that due to the hybridization of the levels near BtransiUon- The 
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increase in AB for T = 10 K compared to T = 0.1 K arises because the transition takes place at lower values of the 
applied field at T = IQ-^ K, ~ 2.60 T, instead of ~ 2.68 T at T = 0.1 K, and because the effect of the Dij vectors 
is larger at lower values of the applied field once the energy levels are close. The matrix elements associated with 
the Dij vectors are larger compared to the Zceman energies in the effective Hamiltonian. As a result increasing the 
magnitude of Dij vectors leads to a larger effect on AB at lower temperatures. It was also observed that the results 
do not depend on the sign of the DM terms. 

The above results were calciilatcd with a gyromagnetic ratio g = 2.0023. A similar analysis for experimentally 
determined temperature dependence of the magnetization at T = 0.1 K gives AB = 0.77 ± 0.10 T and Btransition = 
2.84 ± 0.05 T [18]. In order to compare with this data, an average of the values found from EPR spectra for the 
gyromagnetic ratio is used [17]. They are = gi, = 1.95 and g^ = 1.98, where a, b and c are the crystallographic axes 
of the molecule. Using g = 1.97 and = 70 K, Dv = D'- = 0, it is fomid that AB = 0.75 T and Btransition = 2.81 
T, which is within the uncertainty of the experimental value. Therefore, the inclusion of DM terms accounts for the 
rounding of the magnetization at very low temperatures and provides a description of the magnetic anisotropy of V15 . 
In addition, the gap between the doublets and quartet, ~ 3.7 K, is consistent with the estimates given by Chiorescu 
et al. [18] and Cliaboussant et al. [40] for the effective coupling of the spins in the triangle, induced by interactions 
between these and the spins in the hexagons. Both groups found Jo — 2.4 K which gives a gap | Jq — 3.7 K. We 
also found that the tunnel splitting between the two doublets is Aq = 13 mK, in agreement with the ~ 10~^ K 
estimate [36]. This is an order of magnitude smaller than the estimates given by Chiorescu et al. [28] 0.05 K) or 
Chaboussant et al. (~ 0.1 K). We note that the calculated splittings of the quartets, shown in figure 5, are an order 
of magnitude larger than those of the doublets, shown in figure 4, which may account for the discrepancy between 
experimental estimates and our results. 

In arriving at the form and magnitude of the D^j we have relied on a simple model in which the Dij are non-zero 
only on the strongest bonds and it is not obvious how model dependent our results are in the absence of more detailed 
data. The original estimates by Moriya [31] relied on a supercxchangc model where magnetic order is driven by 
strong short-range correlations due to strongly suppressed hopping between sites. Estimates of the -D^ 's then come 
from considering the mixing of orbitals by angular momentum matrix elements. In V15 interactions between ions on 
hexagon sites are mediated by pairs of oxygen atoms. Depending on whether one of these oxygen atoms is linked to 
an arsenic atom or another vanadium atom, the Heisenberg coupling is either 800 K or 150 K. In these circumstances 
it is not clear how reliable estimates are based on the simple tight-binding approach. The magnitude of the Dij 
vectors found here is comparable to the Heisenberg couplings. This is larger than the estimate based on Moriya's 

original approach which gives l^^jil ~ ^ , where Ag is the shift in the gyromagnetic ratio from that of a free electron. 
However, Katsnelson et al. [39] also found strong Dzyaloshinskii-Moriya interactions with values comparable to the 
isotropic constants in their analysis, which was based on an effective eight spin model of neutron scattering and EPR 
data on Mni2 . 

The determination of Dij is carried out in the absence of the symmetric anisotropy term, Tij. We investigated 
the effect of including Fj j's with a range of magnitudes for its elements with different Dij's. Since the l5y 's are first 
order in the spin-orbit interaction while the Fij's are second order, we have assumed the magnitudes of the symmetric 
anisotropy terms to be smaller than the antisymmetric ones. For illustration we show in fig. 11 the dependence of 
Btransition SLTid the Corresponding values of AB for different forms of Fij given in table 11. In this case D^ = 70 K, 
D^ = D^ = with g = 2.0023. One sees that the values of Btransition and AB are not significantly changed from the 
values given by the DM terms alone. It was observed here that the results depend on the sign of the Fij terms. 

B. Magnetic Relaxation 

Apart from equilibrium magnetization measurements, the magnetization of V15 in a varying external field has been 
measured. The data shows a so-called "butterfly" hysteresis, where the magnetization initially increases with the field, 
then reaches a plateau, and eventually approaches saturation [18,28]. This has been attributed to non-equilibrium 
distribution of phonons, and the spin-phonon mechanism for magnetic relaxation is called "hole-burning". Here the 
role of DM terms in magnetic relaxation is investigated. 

In order to calculate the magnetization in a swept field as in the experiments it is necessary to calculate the 
density matrix for the eight low energy states. However, this is a prohibitively expensive numerical calculation if 
one is to reach applied fields of 7 T with the sweeping rate of 0.14 T/s. This is because the characteristic time. 
To = = 7.64 X 10~^^ s, associated with anisotropy terms is so short compared to the ~ 50 s required to reach 
7 T, that the number of time steps needed is very large. The presence of rapidly varying phases in the elements of 
the density matrix requires a very fine time step <C To- The computational effort to get accurate results is enormous. 
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Therefore, rather than consider a time dependent magnetic field which changes at a constant rate, we consider the 
case in which a 10 T field along the z direction is suddenly switched on, off, or reversed. In the absence of magnetic 
anisotropy there is no change in < 5*^ > with time. However, this is not the case in the presence of DM terms. As 
can be seen in figure 9, there is hybridization between the zero field states as a function of applied field, and matrix 
elements responsible for this lead to a time dependent < > when the applied field is suddenly changed. 

We illustrate this point by calculating the magnetization for the case where the magnetic field has an initial strength 
of 10 T, and = 80 K, = 30 K, = and g = 2.0023. The ground state magnetization is < 5^ >= 1.49. 
Then the magnetic field is suddenly switched off, and the evolution of the magnetization is calculated in the basis of 
the eight lowest states of the same problem in zero magnetic field. Any excited state higher in energy is separated 
by a gap of order J from the lowest energy manifold, and the probability of the state of the system belonging to this 
manifold is 0.9997, so this approximation is well justified. 

The magnetization varies harmonically with time, as seen in figure 12, and the frequencies of its variation are 
determined by the energy differences of the system at zero magnetic field. After a time equal to IOOtq — 10^^° s, the 
field is turned back on to 10 T and < 5*^ > now varies with different frequencies and amplitude. The presence of more 
than one frequency in the time dependence points to the inadequacy of a two-level model and suggests that there 
are significant corrections to any Landau-Zener treatment of the relaxation. The average value of the magnetization 
is different in the two cases, and in particular the average value after the field is turned on again depends on the 
wavefunction at the exact moment the field is switched. The origin of this effect is the magnetic anisotropy of the 
DM terms, a source of magnetic decoherence for the molecule. The calculation can also be done starting from the 
zero field ground states in a non-zero magnetic field, and after ^ 10" i" s switching off the field suddenly. The time 
dependence is similar to the previous case, and once again the average value oi < > after the field has been 
switched off depends on the state of the system at the moment it was switched off. The DM terms therefore provide 
a source of spin decoherence with a decoherence time ^ lO-^^ s. 

In both cases the calculation was done assuming instantaneously switched on and off fields and demonstrates 
magnetic decoherence and hysteresis in the presence of DM terms. If a time varying field is introduced in the 
calculation, the presence of rapidly varying phases requires a very fine time step for the calculation of the time 
evolution of the magnetization. This time step has to be a fraction of lO-^^ s, and a calculation of the evolution of the 
magnetization for realistic times would require an enormous amount of computational effort. However, the small time 
scale calculation demonstrates the intrinsic magnetic decoherence in the presence of magnetic anisotropy. Hysteresis 
is seen in experiment and attributed to the spin-phonon interaction [18]. The present mechanism is independent of 
temperature and can be expected to provide the dominant relaxation process at very low temperatures. 

V. CONCLUSION 

Degenerate perturbation theory was used to study magnetic anisotropy in the molecule V15. The DM term in- 
troduced in the model quantitatively accounts for the rounding of transitions in the magnetization at very low 
temperatures, their width at midpoint, through the avoided level crossing, as well as observed zero field splitting of 
the lowest four states. 

We have shown that the DM interaction is also responsible for spin decoherence and magnetic hysteresis. This is 
present even at very low temperatures where thermal mechanisms are absent. This source of decoherence may pose 
a problem for the application of magnetic molecules in quantum computing, since it is necessary to control the time 
evolution of states with external magnetic fields and electromagnetic pulses [25]. In order for a molecule with a DM 
interaction to be suitable for quantum computing it would be necessary to carry out operations on a time scale much 
shorter than where AE is comparable to the magnitude of the DM interaction. In the present example this is 
10-12 s. 

The use of perturbation theory to calculate the low energy and temperature properties of V15 is an alternative to 
exact diagonalization which has been discussed by Honcckcr et al. [16] for Fcg and Fcg. Perturbation theory works 
particularly well in the case of V15 because of the strength of the hexagonal bonds. However, we have also applied 
perturbation theory to a system of 20 5 = | spins by calculating the effect of quantum fluctuations around the 
classical ground state [24]. The dimension of the Hilbert space in that case is almost exactly the same as Fcs, ^ 10^. 
It is possible that this approach will also succeed in systems with larger, more classical moments and larger Hilbert 
spaces. This would allow a more microscopic approach to the treatment of Mni2, which may resolve some of the 
questions raised by Zhong et al. [44]. 
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TABLE I. Magnetization transition. Index A refers to T = 0.001 K, while B to T = 0.1 K. 
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TABLE IL Magnetization transition with inclusion of Tij terms. = 70 K, = Z)"^ = K, T = 0.001 K. 



point 


ri.i (K) 




ri.2 (K) 




ri.2 (K) 


(K) 


AB (T) 


(T) 


1* 




















0.86 


2.66 


1 














50 





0.69 


2.69 


2 











50 


50 





0.71 


2.69 


3 











50 


50 


50 


0.72 


2.66 


4 








50 


20 


20 


20 


0.71 


2.60 


5 


10 


10 


10 


10 


10 


10 


0.77 


2.59 


6 


10 


10 


10 








20 


0.80 


2.58 


7 











20 


20 


20 


0.79 


2.66 


8 











10 


10 


10 


0.81 


2.67 


9 











10 








0.84 


2.67 


10 











50 








0.86 


2.67 


11 











50 





10 


0.86 


2.67 


12 











40 





20 


0.85 


2.66 


13 

















30 


0.85 


2.65 


14 

















50 


0.86 


2.64 


15 











50 





50 


0.87 


2.64 



9 



14 




FIG. 1. Space configuration of Vib. 
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FIG. 3. Average energy for the lowest and the highest multiplet in the D^, , parameter space: 
(I) D^=D^=0, D^ o: S =^,a: S =^ (II) D"" , D«=D^=0, o: S = ^, A: S = ^ 

(III) D-=L>^ L>^=0, v: S = i, +: S = I (IV) 0^=0^=0' x: S = ^, *: S = ^. 
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FIG. 4. Deviations from the average energy for the two doublets in the , , parameter space: 

(I) o: D'' ^ Dy = 0, (K) (II) □: (K), = = 
(III) o: = (K), = (IV) +: D"" = = (K). 
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FIG. 5. Deviations from the average energy for the quartet in the , , parameter space: 
(I) o: D"" =0^ = 0, (K) (II) □: (K), D'" = = 
(III) o: = (K), =0 (IV) +: D'' = = (K). 



12 




13 




14 



2.9 



2.85 



b 2.75 
1 



2.7 



2.65 



2.6 



2.55 



' 1 ' 

- 


' 1 ' 1 ' 




1 ' 1 ' 1 ' 


'1 


1 

- 


- 

- 











- 


- 


EI 








- 


























, 1 , 


1 1 1 1 , 




1 1 1 1 1 , 


,1 


1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

AB(T) 



1.1 1.2 1.3 



FIG. 10. Midpoint BtranaiUon voisus width AB of the magnetization transition. 

(I) Dy=D'=Q, left to right: = 0, 25, 50, 80, 100 K, o: T = 10"^ K, □: T = 0.1 K 

(II) D'^=100 K, left to right: D^=£)«=0, 20, 30, 40, 50 K, o: T = 10"^ K, x: T = 0.1 K. 
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FIG. 11. Effect of Fij on BtranaiUon and AB. Each point corresponds to a different Fij according to table II. Fjj = for 
point 1*. 
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FIG. 12. < S*^ > vs. time, time in units of tq = = 7.64 x IQ-^^ s, T = 0, = 80 K, = 30 K, = K: "upper" 
solid line: after the field is switched off, "lower" solid line: field switched back on to 10 T. The black circles indicate < > 
at the moments where the field is switched off and back on. 
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